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Abstract 

We compute the reduced cosmic shear up to second order in the gravitational potential without 
relying on the small-angle or thin-lens approximation. This is obtained by solving the Sachs equation 
which describes the deformation of the infinitesimal cross section of a light bundle in the optical 
limit, and maps galaxy intrinsic shapes into their angular images. The calculation is done in the 
Poisson gauge without a specific matter content, including vector and tensor perturbations generated 
at second order and taking account of the inhomogeneities of a fixed redshift source plane. Our final 
result is expressed in terms of spin-2 operators on the sphere and is valid on the full sky. Beside 
the well-known lens- lens and Born corrections that dominate on small angular scales, we find new 
nonlinear couplings. These are a purely general relativistic intrinsic contribution, a coupling between 
the gravitational potential at the source with the lens, couplings between the time delay with the 
lens and between two photon deflections, as well as nonlinear couplings due to the second-order 
vector and tensor components. The inhomogeneity in the redshift of the source induces a coupling 
between the photon redshift with the lens. All these corrections become important on large angular 
scales and should thus be included when computing higher-order observables such as the bispectrum, 
in full or partially full-sky surveys. 
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1 Introduction 



In the early nineties, cosmic shear was predicted to be a promising way to study the distribution of 
matter in the Universe [H [21 13] and since its first detection [H El El [7] it has shown to be a precious 
mean of investigations of the large-scale structure of the Universe, enabling us to explore dark energy 
properties or uncover signatures of mode coupling effects [H El [10] . 

So far cosmic shear surveys have covered only a limited field in the sky. For instance, the 
CFHTLS0 which has produced very promising results over the last years (see e.g. [11] for a recent 
account of these observations), is limited to about a 170 squared degree range. With the demonstra- 
tion of the robustness of cosmic shear observations, (nearly) full-sky surveys such as Pan-STARRS, 
DES, LSST, JDEM, or Euclid are under preparation. They will open the way to new types of stud- 
ies. Akin to CMB observations, such surveys will be an excellent tool to explore the physics of the 
Universe at scales comparable to the Hubble radius, therefore testing genuinely general relativistic 
effects^ In particular, the study of mode couplings, already well established on Newtonian scales, 
can be extended at these very large scales therefore testing the details of our understanding of the 
origin and formation of the large-scale structure. 

Such an investigation requires that we know the types of mode couplings that are expected to 
be seen at such large scales. Calculations have been undertaken to predict the nonlinear growth of 
metric and density fluctuations after modes reenter the Hubble radius |131 [HI [T5] . In the context 
of the CMB anisotropics, progress has been recently made in understanding the effect of these 
nonlinearities, from concentrating on the large angular scales |16 |. ITH] to the details of the physics 
of recombination (see for instance [19} [20] ) . 

So far the investigations of mode couplings in weak lensing were limited to small angular scales, 
corresponding to scales much smaller than the angular diameter distance at the source. Accidentally, 
this distance roughly corresponds to the Hubble radius at the source. Thus, on these scales one can 
consistently neglect general relativistic effects that are suppressed by the ratio between the scale 
probed and the Hubble scale. On small angular scales the dominant contribution to the cosmic shear 
comes from fluctuations of the gravitational potential transverse to the line of sight. Perturbations 
along the line of sight average out and do not yield appreciable effects. In this regime the dominant 
geometrical mode couplings were identified more than a decade ago in |21] , They include the Born 
correction and the lens-lens coupling. In the so-called Born approximation one integrates the lensing 
distortion over an unperturbed photon path. One can consider the correction due to the fact that the 
photon path is perturbed. The lens-lens coupling consists in the correction due to the deformation 
of a distant lens caused by a foreground one. The consequences of these effects have been extensively 
described in the literature and they have been found to have an impact on both the shear power 

1 http://www. cfht.hawaii.edu/Science/CFHTLS 

2 See [TJ] for a recent account of these effects on galaxy clustering observations. 
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spectrum and higher-order statistical observables such as the bispectrum [521 129 IM1 E51 EE E7J . As 
the shape distortion probes the reduced shear rather than the shear itself, there is another correction 
associated to the nonlinear conversion between these two quantities |28[ [29l [27] . Finally, another 
nonlinear effect is the source-lens clustering, due to the fact that the source of the lensed light is 
itself a perturbed field with specific clustering properties correlated with the lens [301 [31]. F° r the 
current surveys restricted to a limited angular field all these types of couplings are undoubtedly the 
dominant ones. 

In view of full-sky surveys one needs to go beyond the small-angle approximation and probe 
scales of the order of the angular diameter distance to the source. In this case fluctuations along the 
line of sight are not negligible and terms other than those described above may become important. 
Furthermore, as we are probing scales comparable to the Hubble size, one needs to undertake a 
full general relativistic treatment. This is important in order to compute accessible higher-order 
observables in full or almost full-sky surveys. In particular, in such surveys, it becomes necessary if 
one wants to compute the lensing bispectrum in the squeezed configuration, when one of the scales 
probed is taken to be much larger than the other two. 

We present here the exhaustive calculation of the weak lensing cosmic shear at second order 
including all general relativistic contributions and without relying on the small-angle or thin-lens 
approximations. However, we do not include in our study the effect of source-lens clustering |30[I31| 
and other intrinsic effects in the alignment and ellipticity of galaxies. We work in the so-called 
generalized Poisson gauge without specifying the matter content of the Universe. We assume that 
there are no primordial vector and tensor perturbations. However, we will take into account vector 
and tensor components of the metric generated at second order from scalar fluctuations. In this gauge 
we derive the reduced - i.e. observable - shear by solving the Sachs equation [32] which describes the 
distortion of the cross section of an infinitesimal bundle of light rays in the geometric-optics limit. 
The advantage of using the Sachs equation instead of the geodesic equation is that it deals only 
with physically observable quantities. As the resolution of the Sachs equation is extremely tedious 
and involves a large number of terms we will develop tests that allow us to check the validity of its 
solution. In particular, some of the contributions to the second-order shear that we compute - and 
that are usually neglected in the small-angle approximation - can be compared to those expected 
from the lensing shear at linear order in a universe with spatial curvature. 

Using the solution of the Sachs equation we will compute the reduced shear by adding the 
nonlinear corrections coming from the relation between this quantity and the shear itself. Finally, 
as in the Poisson gauge hypersurfaces of constant redshift are inhomogeneous, we take into account 
the corrections due to the inhomogeneity of the redshift of the source |33| . The final observable 
reduced shear field that we obtain is a gauge invariant quantity, although its separate contributions 
are not necessarily so. As it is customary for CMB polarization, we express the reduced shear in 
terms of spin-2 operators on the sky. In particular, angular gradients on the sky will be written in 
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terms of spin raising and lowering operators, whose eigenfunctions are the well-known spin-weighted 
spherical harmonics. 

The plan of the paper is the following. In Sec. [2] we give the outline of our calculation. In 
particular, we describe the Sachs equation and how the transverse size of a propagating beam can 
be related to the observable reduced shear including the effects of the source redshift inhomogeneities. 
We will solve the Sachs equation at first order in Sec. while the full second-order calculation will 
be presented in Sec. [H We discuss and comment on our results in Sec. [5] in the context of cosmic 
shear surveys regarding the generation of B modes and the expected contributions to the cosmic 
shear bispectrum. 



2 The weak lensing equations 



We are interested in studying the propagation of a light bundle, i.e. a collection of nearby light rays 
[3H [Ml ESI ESI E3 [38]. We consider two nearby null geodesies a" (A) and a^(A) + £ M (A) that lie in 
the past-light cone of an observer O, connected by a deviation vector field The affine parameter 
A is chosen in such a way that it assumes the same value at O for all geodesies, i.e. Ao = 0. Thus, 
at the observer £^(0) = 0. We denote by k a = dx Q /dA the wave vector of the photons. This obeys 
the geodesic equation, 

Dk^ 

where -Mr = k a V a is the covariant derivative along the geodesic. For rays with an infinitesimal 
separation the connecting vector £^ is also infinitesimal and lies on the null surface, i.e. ^kn = 
everywhere along the geodesies. The evolution equation of the connecting vector reads D^/DX = 
^ u ~V u k^. Indeed, in some arbitrary coordinate system one has 

= k»(x + Sx) - k»(x) + T^k a f 

= ^d^ + T^k a f. (2) 

By taking the covariant derivative of this equation along the photon geodesic and using the geodesic 
equation ([I]) one obtains the Sachs equation [32J, 

^ = Ka^k* , (3) 

where we have used the Ricci identity (V Q V/3 — V^V a )fc' 1 = R^ ua pk u and ~R? VOL n is the Riemann 
tensor. This equation describes the evolution of a light bundle along the geodesic. 

Let us consider the case of a light beam emitted by a galaxy at spacetime position S and 
received by an observer at O. We denote by Vq the observer 4- velocity. It is convenient to define an 
orthonormal spacelike basis nj 1 , with a = 1,2, orthogonal to k^ and to the observer velocity Vq and 
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such that g^n^n^" = 5 ab ^ At the observer position these vectors form a basis, {n^, n^ 1 , k^, Vq}, 
which can be parallel transported along the geodesic, 



DnJ 1 Dv, 



DX °> -D? = °- W 

The subspace defined by {n 1 M (A), n 2 M (A)} is called the screen adapted to Vq and k^. 

We can write the deviation vector £^(A) in this basis. Using the fact that ^k^ = along the 
geodesic, the component along Vq vanishes and we have, for all A, 

with £ a (0) = and f°(0) = at the observer. We can then plug this decomposition into the 
Sachs equation Q. Using the symmetry properties of the Riemann tensor the right-hand side of 
this equation becomes R fJ- ua ^£ i a nJ 3 k u k a . Furthermore, using the equation of parallel transport for 
ria 1 and Vq, eq. flU), and projecting the Sachs equation along the spatial basis nj 1 , one obtains an 
evolution equation for £ a , 

= K\e , (6) 

where the 2D tensor lZ ab is defined by 

n ab = R^pvkrtfnfnf . (7) 
As eq. (|SJ) is linear with initial condition £ a (0) =0, its solution can be written in the form 

e = v ah e b , (8) 

where 



(9) 

A=0 

is the (vectorial) angle between the photon geodesic and the neighboring one at the observer. D ab 
is a linear matrix, called Jacobi mapping. It relates the angle of observation h o to the image on 
the screen adapted to Vq and k^ described by the two spatial components of £ a . From eq. (|HJ), 
eq. ([6]) can be rewritten as an evolution equation for V ab [351 E3 EH [3H] > 

d 2 

j^V ab = TZ ac V cb . (10) 

Equation © and £°(0) = imply that V ab (0) = and d£> afe /dA| A=0 = 5 ab . Note that K ab is 
symmetric but T> ab is generally not. 



3 More generally, it is possible to introduce an induced 2D metric on the subspace described by njf by 
6 



imposing g^vn^n^ = 2T> g a b- This is particularly convenient when employing spherical polar coordinates to 



describe the 3D space. In this case the Latin indices a,b, . . . are raised and lowered by the metric 2D g a b- As 
here we take 2T> g a b = S a b, upper or lower positions of the indices are irrelevant and repeated indices represent 
a summation over a = 1,2. 
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We can decompose the linear mapping T> a b into a spin-0 component and a spin-2 component, 
respectively, 



V = V u +V 22 + i{V 12 -V 21 ) , (11) 
2 V = V u -V 22 + i (©12+D21) . (12) 

The spin-2 part, 2 T> , is directly related to the usual shear spin-2 field 7 = 71 + 172. Indeed, by 
defining the vectorial angular position at the source as 9g = t; a /\s one finds the following relations 

m I I 

7i = -^(V u -V 22 ) , 72 = - — (V 12 + V 21 ) . (13) 

The spin-0 part, qT>, contains a real part, the trace T> = T>\\ + T> 22 , which is related to the usual 
convergence by 

« = l-^D. (14) 

The imaginary part comes from the fact that, unlike lZ a b, T> a j, is not necessarily symmetric. Indeed, 
as we will see, it is not symmetric at second order. The imaginary part of qT> corresponds to a 
rotation of the observed object and it is related to the usual rotation parameter u by 

w = -7jr-(X>i2-X>ai). (15) 
IAS 

At second order, the rotation has no observational consequences on the observed galaxy polarization^ 
Therefore we will ignore it in the following. We are thus left with one scalar degree of freedom 
describing the convergence and 2 degrees of freedom for the shear. The latter can be mapped into 
the so-called "electric" and "magnetic" modes (see Appendix [B]) . Note that at second order the 
electric mode is not necessarily equal to the convergence field as in the linear case. 

What we observe is the ratio between the anisotropic and isotropic deformations, i.e. the reduced 
shear, defined as [39] 

«-T^- (16) 

From eqs. (fT3j) and (fTijl this is given, in terms of the Jacobi mapping V a b, by 

g = - 2 V/V . (17) 

In the following we will solve eq. (|10p in a perturbed universe. For the background we will assume 
a flat Friedmann-Lemaitre- Robertson- Walker (FLRW) metric given by ds 2 = a 2 (r])(—dr] 2 + dx 2 ), 



4 More precisely, one can define the observed complex shape polarization as p = (run — 77122 + 2imi2)/(TOii + 
TO22) where rriij is the luminosity distribution matrix of observed galaxies. For an unpolarized source and in 
absence of rotation this is p — 2g/ (1 + gg*) [39] . In case of rotation this relation becomes (1 — R*)/(l — k)p = 
2g/(l+gg*) where K is complex and defined as k = k + iuj and g = g/(l — k). Thus, as k is real at linear order 
and 7 vanishing at zeroth order, the imaginary part of k enters in the expression of the observed polarization 
at third order only. 
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where r\ is the conformal time. Since null geodesies are not affected by conformal transformations, it 
will be convenient to perform the calculation without the Friedmann expansion and reintroduce the 
effect of the expansion only at the end. Indeed, as shown in Appendix^] the effect of the expansion 
can be simply taken into account by rescaling the mapping T> a b by the scale factor, e.g. T> a f, — > aD a b- 
Note that, as both 2^ and T> get rescaled by the conformal transformation, the reduced shear (|17|) 
is not affected by the expansion. We parametrize the photon geodesic such that x°(X(r])) = r/o — 77. 
It is thus convenient to define x = Vo — V so that 



and the evolution equation ([10]) now reads 

d 2 1 dA; d 1 , , 

d^ + ^d^d^^ = (fcop 7 ^ 2 ^ • ( 19 ) 

We can solve this second-order differential equation order by order in the metric perturbations. 

The solution of eq. (|19p gives the linear mapping between the observed angle and the shape of 
a source at a given coordinate time rjs- Observationally, we are interested in a mapping where the 
source is defined at a given redshift z$. As in an inhomogeneous universe the redshift is a perturbed 
quantity, at second order we expect a contribution to the reduced shear due to the coupling between 
this perturbation and the lenses. Thus, the reduced shear at constant redshift is given by 

g z = g--^-Sz s , (20) 
dz s 

where 5z = z + 1 — a^/a is the perturbation of the redshift. As this correction is not conformally 
invariant, it introduces a dependence on the expansion. At first order, though, only the trace of 
the linear mapping, V, is affected by redshift perturbations [33], but the traceless part is not and 
therefore in this case g z = g. 



3 The shear and the convergence at first order 

As a warm up exercise, before the calculation of the shear at second order we derive here, using 
eq. (|19p . the shear and the convergence at first order. We consider a perturbed FLRW metric in 
Newtonian gauge, written in Cartesian coordinates as 

ds 2 = a 2 (r]) [-(1 + 20)dr/ 2 + (l - 2i/>)dx 2 ] , (21) 

where we have neglected primordial vector and tensor perturbations. For convenience, we define 
also the Weyl potential VP as 

tf=(0 + ^)/2, (22) 
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and we will use it in the following whenever the combination + ip appears. As explained in the 
previous section, the reduced shear (fT7|) is not affected by the expansion and we can set the scale 
factor a = 1. 

Let us define £;^(0) = (l,e r ) as the photon 4- momentum at the observer, where e r defines the 
direction of the line of sight. Note that in this definition we have set the metric perturbations at the 
observer position to zero. The final result will be independent of this choice. Indeed, since metric 
perturbations at the observer do not depend on the direction of observation, they can be reabsorbed 
into the homogeneous mapping. As in flat spacetime the Christoffell symbols vanish, k^(0) is parallel 
propagated along the background geodesic, while the curvature tensor R^afiv vanishes thus making 
lZ a b at least a first-order quantity. Then, at first order in perturbations eq. (|19p simplifies to 

3— T^T> a b + — — -7—'D a b = 1Zac'Dcb, (23) 

a-X 2 dx dx 

where we have used that dk° /dx vanishes on the background. 

On the background this equation becomes d 2 V ab /dx 2 = 0. Requiring that the homogeneous 
mapping is proportional to the unit matrix yields T> ab = x$ab for the background solution. Further- 
more, we can plug the background solution for T> ab in the second term on the left-hand side and on 
the right-hand side of eq. ([23j) to obtain 

d 2 ^ dk° 

dx 2 dx 

We then integrate this equation once up to the source using the boundary condition dV ab /dx\o = $ab, 



172 V ab + — d ab = XK ab ■ (24) 



— V ab = (2 - k°)5 ab + / d X 'x' K ah . (25) 
dx Jo 

The solution of this equation, after an integration by parts, can be written as 

r\s rx.s 
V ab {xs) = / dx(2 - k°)5 ab + / d X (x5 - x)xK ab ; (26) 
Jo Jo 

k° can be obtained by solving the photon geodesic (P) at first order, 

k° = 1 - 2<P + / dx'2V , (27) 



where the dot denotes a partial derivative with respect to x, i.e. '= d/dx- 

Let us now compute lZ ab . We denote by n a ^(0) = (0, e a ) the spatial basis at the observer. As for 
the 4-momentum of the photon, 71^(0) is parallel propagated along the background geodesic. Thus, 
according to the definition ([7]) , in order to compute lZ ab we need to contract the curvature tensor at 
first order with the unperturbed £^(0) and n a /i (0), which yields 



Kab = -e a V2* y - P^Sab • (28) 
dx 



S 



Plugging lZ a i, given in this equation and the expression for k° given in eq. ([27]) into eq. ([25]) and 
integrating by parts d 2 ip/dx 2 we obtain the Jacobi mapping at linear order, 

rxs r -i . rxs 

V ab ( X s) = (l-4>(xs))xsKb+ / d X 4$- 2( X s ~ X)* <U-e a V / d X (x5"X)x 2*,y . (29) 

Jo L J JO 

According to the Born approximation we can evaluate the integral along the background geodesic 

so that in this expression ^ = ^(x, e r x)- 

Note that the metric (|2ip is conformal to ds 2 = — (l+4^ , )dry 2 +dx 2 . Thus, as photon geodesies are 
conformally invariant, one would naively expect T> a b to depend only on the combination \& = (<f)+ip)/2 
[37] . This is too quick a conclusion. Indeed, parallel transport of the basis nj 1 is not conformally 
invariant and the basis is deformed by the spatial curvature at the source position. For this reason 
the first term of eq. ([29"]) also depends on the curvature potential at the source, ips- 

We now want to extract the shear and convergence from the mapping T> a i,{xs)- I n order to do 
so, it is useful to introduce spin operators on the sphere (see Appendix [B]) . To each point of a 2D 
Riemannian manifold described by an orthonormal basis {ei,e2} we can associate a spin-s field S X 
such that under the rotation of e a by an angle a it transforms as S X —> e ias s X (for more details see 
[40 1 14 H I4"2]). More precisely, following [4T], the local freedom in the choice of the basis is equivalent 
to the transformations 

e± = ei ± ie 2 -»• e iQ e± . (30) 

To every spin-s S X we can associate a symmetric and trace-free tensor of rank s > 0, X ai ___ ag : for 
s > 0, 

X ai - a ° = 2- s s Xe a 2 ■■■e a J . (31) 

The inverse relation is 

s X = e^---e a + s X ai ,„ as . (32) 

For s < we define X ai "- a M = 2-l s l s Xe+ 1 • • • e^ s| . 

As we are interested in describing the lensing field on the sphere of the sky, for the orthonormal 
basis we choose the two coordinate basis vectors of a spherical polar coordinate system, {ei,e2} = 
{e^e^,}. Let us rewrite the three spatial vectors e r and e a in the Cartesian representation as 

e r = (sin 8 cos ip, sin 6 sin 93, cos 9) , (33) 
eg = (cos 8 cos if, cos 6 sin 93, — sin 6) , (34) 
e v = (— simp, cos if, 0) , (35) 

where 8, cp are the angles of observation. Furthermore, we can define operators that increase or 
decrease the index of the spin by 1, 

S X = - sin s 6{d e + icsc 0<9 v )(sin^ 6) S X , S X = - sin~ s 8{d e - icsc^)(sin s 6) S X . (36) 
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With the definitions ([3"3"1)-([3"5T) above we have 



e\di = d r , e l e di = -d g , e l d { = — — -cL , (37) 
X v xsm6» 

where we have used that r = x along the photon geodesic. Using these relations it is easy to verify 
that, if X = X(x) is a scalar field, 

eiXi = --0X, elXi = --$X . (38) 
X X 

Then, employing the useful relation x e ±d% e ± = c °t one can verify that 

< .'< ^X,j = p 2 X , ejeiXy = ^fx, (39) 
and, analogously with x e ^i^± = —cot9e± — 2e r J , that 

e^eiXy = eie{X yij = \ 0X - -X >r . (40) 

X X 

Note that the 2D Laplacian on the sphere is given by 0X = 0X. 

Let us first apply these definitions to TZ a b- By making use of eq. (f32|) we can define a spin-2 field 
on the sphere, 2^ = e+e+lZ a b, which by eq. ([39]) reads 

2 

X 

Then, using that e^ed = (e+ej? + elel)/2 and eq. (jiT)]) . the trace of TZ a b reads 



2 n = -^0H. (41) 



K = -U,.-- 2 00*-2p^. (42) 
X X 2 0.x 2 



We can do the same for the mapping (|29p . One finds the spin-2 mapping field, 

2 V( XS ) = -2 r d X ^^ 2 * • (43) 
Jo X 

At first order the reduced shear g is given by 2D divided by the background part of T>, which is 
simply V = 2xs, i-e. 

g = [ XS d X ^^ 2 ^ • (44) 
Jo XXs 

Taking the trace of eq. (|29p and defining 5T> = T> — 2x one finds, after an integration by part, 

8V{ Xs ) = -2iP{xs)xs + 2 £ & 'dxU* -M^A00^\ ) (45) 

which is proportional to the convergence, K. The first two terms on the right-hand side are not usually 
included in the convergence because they are negligible on small angular scales f] In particular, the 



5 The Laplacian $ $ introduces a factor 1(1 + 1) in harmonic space, thus enhancing the contribution to the 
power spectrum from the third term of eq. (|45[) with respect to the contribution of the other terms. 
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first term is a relativistic effect due to the deformation of the size of the source caused by the 
curvature potential at the source. The second term is the Shapiro time delay [13]. Both effects are 
negligible on small angular scales. 

We can check that the first two terms in eq. (|45p contribute to the convergence by considering a 
nonexpanding FLRW universe with nonzero constant spatial curvature K and metric 



1 _, . o i ,„o Kr 2 

,2 ( 



ds 2 = -dri 2 + - ^(dr 2 + r 2 dQ 2 ) with ip K (r) = ^— . (46) 

(i+#W) 2 



In this universe the angular diameter distance is 



DM = [**SK*ISK fa A->0, 

1 sinh(^i^x)/^^l for K < . 



For small curvature, K <C 1, the expression above can be expanded at first order in K to give 

Da(xs) = XS-^L. (48) 

The first term on the right-hand side is the angular diameter distance in a spatially flat universe 
while the second term is a small perturbation to it, due to the curvature. 

As the trace T> corresponds to twice the angular diameter distance, we can check that 5T>(xs) /2, 
with 5T>(xs) given in eq. ()45|) . correctly reproduces this small perturbation in this particular case. 
For = and ip = ipx (and a = 1) the flat perturbed metric in Poisson gauge, eq. (f2T|) . reproduces 
at first order in tpK the spatially curved metric (|46p . Thus, in this case * = (0 + V)/2 = Kr 2 /8. As 
the Weyl potential ^ depends only on the radial coordinate, the last term on the right-hand side of 
eq. (|45p vanishes so that it does not contribute to the angular diameter distance. Furthermore, the 
contributions from the first and the second term of eq. (|45p are, respectively, — K\%/^ and 
where we have used that, at lowest order in vpx, r = X along the photon geodesic. This yields 
ST>(xs)/2 = — i-e. the second term of eq. (@8|). Note that the first two terms of eq. ([43]) 

have no counterpart in the shear. Indeed, a constant curvature deviates light rays only isotropically. 
Furthermore, Kaiser's relation [44] between shear and convergence, 0k = 0^/, is only valid in the 
limit of small angles, when the first two terms of eq. (|45p are negligible. 



4 The shear at second order 



In this section we will compute the lensing shear at second order. As we concentrate on the shear, 
we will compute only the trace- free and symmetric part of the matrix T> a b- Indeed, at second 
order the trace and the antisymmetric part contribute only to the convergence n and the rotation 
u, respectively. We will do the calculation in the so-called generalized Poisson gauge, where the 
second-order metric can be written as 



ds 2 = a 2 ( V ) -e^dr? 2 + 2w< dr/dx i + (e" 2 % + hij) dx i dx j 



(49) 
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Here the vector component cjj is divergenceless, diUJi = 0, and the tensor component hij is diver- 
genceless and traceless, dihij = = ha. As we neglect primordial vector and tensor perturba- 
tions, u)i and are only second-order quantities. Note that we have used the exponential form 
for the gravitational potentials in the metric. Indeed, in this form the metric is conformal to 
ds 2 = — e^drj 2 + 2wj drjdx 1 + (5{j + hij) dxMx- 7 so that the effect of scalar perturbations on the null 
geodesic is only through the Weyl potential \P = {<p + ip)/2. As done in the previous section, we will 
drop the effect of the expansion setting a = 1. We will reintroduce the expansion in Sec. 14.31 



4.1 Solving the Sachs equation 

To solve the Sachs equation at second order let us go back to eq. (fl"9|) . This can be integrated to 
obtain 

ic Vab = 6ab + £ dx ' {vh nacVcb ~ ^k Vab ) ' (50) 

where we have used that dV a b/dx\o = $ab at the observer. Since dk® /d\ vanishes on the background 
we can rewrite the last term of this equation using the first-order equation (|25p and integrate eq. (|50p 
up to the source S. We obtain, after integration by parts, 

VabiXs) = XS^ab + / ^ S ab , (51) 

Jo X 

where we have defined the source term S a b as 

S ab = j^n ac V cb - |^(2 - k°)5 a b - x^- £ dx'x'Kab • (52) 

At second order we need to go beyond the standard Born approximation. Since the source term 
Sat, is at least first order, we just need to integrate along a geodesic which is perturbed at first order. 
Thus, we evaluate the source term at the perturbed geodesic position x l peTt (x) = x l {x) + ^ x% {x)i 
where 5x' 1 is the geodesic deviation. Expanding S a b along the background geodesic, 

Sab{x l pcrt ) = <Sab(x l ) + Sx J • 5(S a b)j\x , (53) 

where the shift 5(S a b)j will be computed below, we can rewrite the source term S a b as a function of 
the unperturbed geodesic position x 1 as 

X 2 X dk° dk° f x 

Sab = 7^0j2^ a6 + X^ac SV cb ~ ^0 "^( 2 ~ k °)$ab ~ X~^ J dx'x'Rab + X' ■ 5{lZ a b)j , (54) 

where to write the first two terms we have separated T> a b into its background and first-order part as 
T^ab = X$ab + ^ab and we have employed for the last term that at leading order S a b = x 2 ^ab- We 
will now compute the symmetric traceless part of the source S a b evaluating one by one each of the 
terms on the right-hand side of eq. (|54p . 

Let us start with the first term in eq. (|54p . Since on a flat background spacetime the unperturbed 
part of R^afiv vanishes, according to the definition 0, i.e., !Z a b = R^ai3uk a k^ nj 1 n b u , in order to 
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compute lZ a b up to second order we just need to consider the photon wave vector and the basis 
ria up to first order. Integrating the geodesic equation (P) one obtains, up to first order, 

k° = l-2(f>+ dx' 2^ , (55) 
Jo 

f'X 

= (1 + 2ifj) e ; - / d X ' 2d 1 ^ . (56) 
Jo 

Integrating up to first order the parallel transport equation for nj 1 , eq. (HD, one finds 

n a ° = -e* [ X d X 'd t 4>, (57) 
J o 

n* = efr + ejej ' • (58) 
Jo 

Since the vector and tensor parts of the metric are already second order, they enter linearly in R^a/Su, 
which just needs to be contracted with the unperturbed k^ 1 and nj 1 . Furthermore, it is convenient to 
compute the ratio k l /k° which appears when combining k l in the definition (J7J) with the denominator 
of the first term of eq. (|54p . This yields 

= e/(l + 2*) - ejej J dx' 2*,, . (59) 

Finally, computing R m p v up to second order and making use of the definition (J7|), the traceless 
part of the first term of eq. ([54"]) reads 

X, {n ab ) T = - x 2 efaj f 2* - + 4*^- + 8** - + 8e r ** « |* 

j . , 2 . °, N (60) 

k _d_ 1 fc . 1 r , l_d_ fc \ 

— e r + ~^ UJi >i ~~ 6r kl M ~ 2 dx 2 " 2 dx^' feij i ' 

where the parenthesis in the indices denote symmetrization. After few manipulations, employing 
eqs. (|38p and (|39p . and the useful relation x e +d% e r = e+, we can write its contribution to iS = 

- 2 2 tf - 4 0tf 0* - 8* <f ^ - 8 (x 0V, r ~ 0V) / dx'— 0V 



e?e}S nh as 



+ f^r + ^(X l") + ~ # 2 /w + ^-(x iM + | ^ (X 2*0 



2 ^ (61) 



where cj r = e r l Ui, h rr = e r l e r hij and we have defined the spin-1 part of lj, \oj = eluji, the spin-1 
quantity %h r = e+ej and the spin-2 part of the tensor mode hij, 2^ = e+e^/ijjlfl 



6 Note that uii and hij can be decomposed respectively as = u) r e r l + e_|? + |iw ej 8 and = 

/i rr {e^ej — ^e^e.) + _i/i r e+ef^ + \h r e_e^ + \-ih z+ e + + 52^ eje_£. These components are not 
independent, as they are related by the divergenceless conditions u>ij = = /i^. 
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The second term of eq. ([53]) is the product of two first-order quantities, lZ ac and 5V cb , both 
symmetric at first order and containing a trace and a traceless part. The symmetric traceless part 
of this product is simply given by 

X (n ac 5V cb ) ST = X (K ab ) Tl -5V + X \n{5V ab ) T , (62) 

where 5V is the perturbed trace of the mapping V ab and 1Z is the trace of 1Z a b- Note that the term 
X^a C SDrf, contains a trace and an antisymmetric part and thus contributes only to the second-order 
convergence and rotation. Using the expressions for 2TZ and 1Z, eqs. ([31]) and ([32"]) . and those for 
and V, eqs. ([35]) and ([35]) . one obtains for the contribution of the second term to 2*?, 



(63) 



2V> 2 ^ - - 2 ^ f X dx' (2* - - — — 0^ 
X Jo V X 

+ 2(m tr + ±.009 + xp£) f X &x'^^ 2 * 
V X dx V y X 

The first line is x while the second is xR- 

The third term of S ab is a pure trace so that it does not contribute to the shear. The traceless 
part of the fourth term of eq. ([53]) can be straightforwardly computed by noting that 

dk° d0 ■ 

— = -2-f- + 2tf . 64 

dx dx 

Thus, its contribution to the source 2^5 is 

4x(-^ + ^)^dx'^^. (65) 

In order to express the last term of eq. (|54p in terms of a spin-2 field we need to solve the geodesic 
equation at first order. We can solve dx l /dx = k l /k° making use of eq. ()59|) . After integrating by 
parts this yields the geodesic deviation, 

5x l = e/ f X d X ' 2* - ejej ^ dx'ix ~ x') . (66) 
Jo Jo 

The shift of !Z a b, 5{TZ a b)i, is simply given by the variation of \P and d 2 V>/dx 2 along the geodesic in 
eq. ([28]) . We only need to take the traceless part of 5(JZ ab )i. Thus we have 

S(n ab )J = -eje b k , (67) 

which can be contracted with 5x l of eq. (j66p . Note that, as we are varying directly the scalar Weyl 
potential \E f , it was not necessary to introduce a covariant derivative on the sphere as in [46] , With 
the definitions and relations of Sec. [3] it is possible to verify that 

ele{e*X tijk = ~ 3 X , e^e^X >ijk = -±^00*X-^(^ 0X, r ~^0* 
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Using these relations and e^e a = (e+e_ + e^e+)/2 one finds, for this last contribution, 



4 - 



d X '* 



+ 



Hp 



0^ + - 00 2 ^ I ' d X '^-^- $^ + 8 ( 0*, 



2 T7 



X' 



X 



x - X 



X' 



X 



0V 



x - x' 



X 

(69) 



where the first and second lines come from contracting with the first and second terms on the right- 
hand side of eq. (|66p . respectively. Using the commutation rule for the spin raising and lowering 
operators [41] . 

(00- 00) s X = 2s s X, (70) 
it is convenient to rewrite the second term in the second line of eq. (|69p as 



X 



- 2 ^ I d x '^^0y = -000^ I d X '^^0y + z 0y I dx'^-^0*. 



X' 



X 



x - X 



-x - X 



(71) 

X' X Jo " X 

Finally, combining eqs. (f6"Tj) . ([63]) . ([65]) . (|69|) and ([7T]) . replacing partial derivatives with respect 
to r by using d r = d/dx — d x (we remind the reader that ^ = d x ^) and integrating the total 
derivatives by parts, we obtain 



iD{Xt 



dx 



XS~X 



dx 



o 



X 

XS ~ X 
X 



0^ + 



1 

X 



dx'^—£- 0^+00-9 I dx'^-^- 0^> 
X ' ' 



1 + l^ (xs) ^ _ ^ I T ^ 02y + ^2 

2 2 x Jo 



/ X X 

X' 



dx'* 



o 



' i% dx'^^0^ 

X 



X Jo 

+ ^2h(xs) + 



rxs 




-4 / d X 


V 


Jo 


Jo 



X 



x.s 



^X 



Xs - X . , 1 L , Xs 



X 



u r + -h rr ) H 0(icj + i/i r ) 

2 X 



(72) 



This is the solution of the Sachs equation. It is written in terms of spin raising and lowering operators 
and 0, which are just the extensions on the sphere of the usual angular gradients: for small angles 
one can replace them by angular gradients V er transverse to the line of sight. The advantage of using 
this representation is that the eigenfunctions of these operators are simply spin-weighted spherical 
harmonics (see Appendix B). 

In the small-angle approximation, i.e. when transverse scales are much smaller than radial scales, 
eq. (|72p reduces to its first line. Indeed, in harmonics space the presence of the operators and 
is associated to an I factor. Thus, when observations are confined to large / one gets a larger 
contribution from those terms containing more and operators, such as the first line. The 
parentheses of the first line contain the usual couplings considered in the literature [21\ [22| [23l [24"1 [25] , 
i.e. the lens-lens correction and the correction to the Born approximation. Note that there are other 
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couplings at play: the nonlinear growth of dark matter fluctuations induced by gravity on small 
scales introduces other second-order effects that are incorporated in the Weyl potential ^jj) 

In the small-angle approximation the last three lines of eq. (|72p are suppressed by the ratio 
between the transverse scales probed and the longitudinal distances. However, for a full-sky survey 
where one probes larger scales they can become important. Let us list here these terms collecting 
them by their physical interpretation (for simplicity, integration over dx and dx' will be omitted): 

• An intrinsic contribution, integrated only once along the line of sight, 

2 XS ~ X 2 V 2 {x) , (73) 
X 

coming from the Riemann tensor at second order. This is a purely general relativistic effect 
of second-order gravity. 

• A source-lens coupling, 

2*^*^ xs ) $H\x) , (74) 
X 

which comes from the coupling between the lens and the curvature at the source, inducing a 
deformation of its shape. As we will see, this term is absent from the final expression of the 
reduced shear, as it cancels with an equivalent term coming from the corrections due to the 
denominator of eq. (|17|) . 

• Time delay-lens couplings, 

A Xs_zX fy^L 02^) + 0H( X )*(x') + *(x) 2 ^(X", 

V x x x J 

which come from the coupling between the lens and the time delay that occurs during the 
longitudinal photon path. 

• Deflection- deflection couplings, 

4 xs_zx L W) + ^ {x) x^X_ _ 8 ^ {x) ^ {x/) (76) 

XV X J X 

which are due to the couplings between two changes in the photon directions. Note that 
the photon deflection is described by a spin-1 field (the deflecting angle), i.e. a spin raising 
operator acting on a scalar. Taken alone it does not change the shear as it just affects all 
photons of the beam in the same way. However, the coupling of two deflections generate a 
spin-2 field which contributes to the shear. These corrections, as well as the time delay-lens 
ones, are integrated twice along the line of sight. 



7 To a large extent, these effects are those that are expected to dominate in the current surveys and that 
lead to detectable effects [2i~ l |47 l |48]. 
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The last line of eq. (|72p contains the effects induced by vector and tensor modes generated at 
second order. As vectors and tensors are already second-order quantities, they enter linearly in 
this expression. The vector component enters only through terms integrated along the line of sight. 
As the tensor component is a spin-2 field it induces also a boundary term, which accounts for the 
distortion of the shape at the source. The integrated contributions from the tensor component agree 
with what was found in |491 150] . Note that, at second order, the separation into scalars, vectors and 
tensors done here is gauge dependent. Indeed, we expect all these different contributions to give 
comparable effects to second-order observables such as the bispectrum. This is similar, for instance, 
to what happens when one computes the CMB bispectrum on large angular scales [14j. Finally, note 
that our final result (|72p cannot be written as the action of 2 on a scalar quantity. Thus, the shear 
will contain also B modes (see Appendix iBj) . 

4.2 Testing the solution 

As the reader has certainly realized, the derivation of eq. (|72p is extremely tedious and involves 
many steps. Thus, it is important to develop tests in order to check this equation and, in particular, 
the new nonlinear couplings (|73p -(|76 p . One implicit check is that 2D behaves as a spin-2 field 
under rotation of the screen basis. This is automatically ensured by the use of the operators 
and instead of the angular derivatives. Furthermore, we can check the terms in eqs. (|73|) -(|75 p by 
studying specific cases where part of our calculation can be pursued nonperturbatively. For instance, 
we can study the shear at linear order in a universe with constant curvature. In the limit where 
the curvature is small we must recover the couplings between the curvature and the gravitational 
potential given in eqs. (|73p -(|75 p . As we will see, such a strategy can be generalized to other cases. 
Unfortunately we were not able to develop an analogous test for the terms appearing in eq. (|75|) . 

4.2.1 Couplings from perturbing the spatial curvature 

Let us consider a perturbed nonexpanding FLRW spacetime with constant curvature K and metric 

ds 2 = - (1 + 2cf) dr/ 2 + 2 ( dr2 + r 2 dfl 2 ) , (77) 

(l + i/j k ) 

with 

cf> = 4>(x») , ^ K = ^ , (78) 

where we have perturbed only the 00 part of the metric. At linear order, the spin-2 mapping field 
2~D for this metric is (see for instance |21j ) 

2V K ( XS ) = - I" d X D ^ S ~ x) ?*, (79) 
Jo D A(X) 
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where Da{x) is the angular diameter distance given by eq. (|4"7|) . Equation ([75]) is just the gener- 
alization of eq. (03]), where Da(x) = X> to a constant curved FLRW universe. Since ipx does not 
depend on the angles, 2 ipK = and thus only 2 (f) appears on the right-hand side of this equation. 

For a small curvature K we can expand eq. (|79p at first order in iftx. The angular diameter 
distance ([37]) reads Da(x) — X ~~ ^X 3 /6- Furthermore, we have to evaluate <fi on the geodesic 
solution for a curved universe, i.e. at r(x) = X + -^X 3 /12, which yields 

PttrMA <P) = 2 HX, 0, <t>) + 2 Mx, 0, 0) . (80) 

Plugging these expressions into eq. (|79|) . replacing the derivative with respect to r using d r = 
d/dx — d x and integrating by parts we obtain, up to first order in ipx, 



••xs 

2 D K = - I d X 

o 



2 <P - % ( ^ - 3x 2 - 6x1 + Qxxs) 2 4> -^(xs- x)x 2 2 4> 

X 12 V X J 12 



(81) 

The last two terms on the right-hand side of this equation can be seen as "second-order" corrections 
to the first-order expression (|43|) . of order ~ 0(4>iI>k), due to the coupling between the gravitational 
potential (j) and the curvature perturbation tpx- These corrections are already incorporated in 
our second-order expression ([72]) . Indeed, by replacing \& = (</> + tpx)/ 2 in this equation, neglecting 
second-order terms of order ~ 0((p 2 ) but keeping those of order ~ O(cj)ipK) one finds, after integrating 
by parts, eq. (f8Tj) , 

This calculation can be extended to a spacetime with radial-dependent curvature, i.e. ipK( r ) a 
generic function of r. In this case the spin-2 mapping 2D is given by 

2V K {xs) = - J** dx ^*^ 2 ^ , (82) 

where the Green's function G(xsi x) an d the angular diameter distance Da{x) can be derived from 
using eq. (fl~9|) in a homogeneous universe with spatial curvature ipK- In this case the trace of eq. (fl~9|) 
becomes 

= Urmbd . (83) 
where, at first order in ipx, R = — ipK r / r ~ rri implying 



1>K = -J dx'/x'j d X "x"R(x"). (84) 

The two solutions of eq. ([83]) . Da(x) an d C(x), are determined through their initial conditions, 
i.e. Da(x) — > X and C(x) -> 1 + 0(x 2 ), so that 

^A(X) = X + [ X d X ' [ X d X " X "R(x") , (85) 
Jo Jo 

C(x) = 1+ [ X d X ' [ X d X "R(x") ■ (86) 
Jo Jo 
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The Green's function is then given by 



G( X S,X) = D A (xs)C( X ) ~ D A ( X )C(xs) • (87) 

Evaluating 2 <p on the geodesic solution for a curved universe, and replacing eqs. (|85p and (JH7J) into 
([82]) we obtain, at first order in <f> and ipx but keeping terms of order ~ Ofyipx), 

/"XS 1 ( /*X /*X 1 /*x /*x' \ 

2V K {xs) = - dx$ 2 4> d x ' d x "R(x")-- d x ' dx"x"R(x")) 

Jo [ X \ Jo Jo X Jo Jo J 

rxs rx' i rxs rx' 

- d x ' d x "R(x") + - d x ' d x "x"R(x" 

J\ JO X Jv Jo 



'x J o A J X 

xs "X.& — X 2 f ^ 

dx 4>,r I ipKdx' ■ 

o X Jo 



We have checked that, after integration by part of the last line, our general expression (|72p reproduces 
this peculiar case. 

4.2.2 Reparametrization invariance under a time shift 

Until now we have tested the couplings between the 00 metric perturbation (j){x^) and a radial- 
dependent spatial curvature ipK(r). Analogously, we can test the nonlinear couplings between the 
spatial metric perturbation %p{x^) and a 00 metric perturbation which depends only on time, 4>t(i])- 
Such a perturbation can be reabsorbed into the time coordinate through a homogeneous shift of the 
time dx = dx(l + 4>t(x))- As a homogeneous time shift does not change the gauge, we expect our 
expression ()72|) to be invariant under a first-order coordinate change 



rx 

X = X + / dx 
Jo 



>T , X = X . 



Let us check that this is the case. At first order in (j)?, under this coordinate transformation 
the metric perturbations change as (ft(x) = 4>{x) ~ 0t(x) an d V'(x) = VKx)- Evaluating on the 
geodesic r = x = X + Jq dx'^T, the first term of eq. (|72|) transforms as 

-2 [ XS d X XS ~ X $ 2 ^> = - 2 [ XS dx XS Z X 2 ^> 



o X Jo X ^ 



rxs i rx rxs -9? — y ■ f x 

2 / d*- 2 * / dx'<f>T + 2 / dx^^ 2 ^ / dx>T • 
Jo X Jo Jo X Jo 



One can check that replacing ^(x) = ^{x)~4>t{x) an d X = X m to the second-order terms of eq. ([72 
exactly cancels the second line of eq. ()90p . leaving invariant under the transformation (|89|) . 
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4.2.3 Reparametrization under conformal transformation 

In Appendix [X] we show that under the conformal transformation of the metric g^ u — > VL 2 (x^)g^ u 
the mapping T> ab transforms as 

v ab ^nv ab . (91) 

If the conformal factor is just f2 = l + <5f2, where <5f2 is a small perturbation of order ^, this conformal 
transformation is equivalent, at first order in 5£l, to a redefinition of the potentials in the metric, 

ip^ip-Sn, + Sn . (92) 

Since, as expected, the Weyl potential ^ = (<fi + ifi)/2 does not change under this transformation, 
in eq. (|72p only the boundary term proportional to tp(xs) is n °t invariant. Transforming this term 
according to (|92p yields a contribution to 2 'D(xs) i n the new metric, 

2 V( XS ) 2 V( XS ) - 2<5fi(4) d X^^ ^ ■ (93) 

Jo X 

This is exactly what we expect from the conformal transformation 2 D — ► (1 + #0)2^ • 

4.3 The reduced shear 

As mentioned in Sec. [2j the quantity that we measure is the reduced shear, which is given by the 
ratio between the spin-2 anisotropic mapping 2 T> and the trace T>, eq. ()17p . Expanding this equation 
up to second order using that V = 2xs + S"D we obtain for the reduced shear 

2 V 2V 5T> 

+ <94> 

where 2 T> in the first term on the right-hand side is given by eq. (|72p and, from eqs. (|43p and (j45|) . 
the second-order correction on the right-hand side is given by 



(2xs) 2 Jo XXS Jo X'XS 

+ ip(xs) r d X ^^ 2 ^ - 2 r d X * r ^'^^ ^ 

Jo XXS Jo Jo XX S 



(95) 



X'Xs 

The first line is the usual correction to the reduced shear due to the coupling between the convergence 
and the shear [28} I29j . The two corrections in the second line are negligible in the small-angle 
approximation. The term proportional to ip(xs) comes from the coupling between the lens and 
the curvature perturbation at the source contained in the trace of T> a b, see eq. (|45p . Note that it 
cancels with the one in the expression of 2 T>. Indeed, as both the isotropic and anisotropic part 
of T) ab change according to eq. (|91|) under conformal transformation, we expect the reduced shear, 
which is their ratio, to depend only on the Weyl potential which is invariant under conformal 
transformation. The last term in the second line is the coupling between the time delay contained 
in V with the lens. These terms are thus of the same order as those discussed in eqs. (fT3|) - (fT5|) 
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The expression of the reduced shear above is given taking the source at constant conformal time 
775. However, in order to relate this quantity to observations we need to compute it using a constant 
redshift z$ for the source, given by 

*8 = ^-l. (96) 

As in the Poisson gauge the z = const hypersurfaces do not coincide with the 77 = const hyper surf aces, 
the redshift zs is not homogeneous and we expect a correction to eq. (|94p coming from the coupling 
between the perturbed redshift plane of the source and the lens. As explained in Sec. [21 eq. ([20]) . 
the reduced shear at constant redshift is given by 

9z = g + 5g z , 5g z = -p?- 6z s , (97) 

dxs dz s 



where 5zs is the perturbation of the redshift (j96|) . We must now reintroduce the expansion of the 
Universe. By taking the derivative with respect to xs of the linear expression for g, eq. (|44[) , we find 

dg dxs 1 f xs j_ 5*2 



d X , (98) 



^XS dz s x s H S 

where H is the Hubble rate defined from the cosmic time t, dt = — adx, as H = . 

To compute 5z$ we can perturb at first order eq. ([96]) using the expression for k° given in eq. ([2j 
Setting to zero the perturbations of the metric and of the velocity at the observer position, we obtain 

5z s = -(l + z s ) U(xs)+e r \ xs -v s -2 £ $ d X ^ ■ (99) 

In the three terms on the right-hand side of this equation one recognizes the Sachs- Wolfe, the 
Doppler and the integrated Sachs- Wolfe effects contributing to the photon redshift perturbation. 
Finally, combining the expression (|98|) and the redshift perturbation (|99|) we obtain for the redshift 
correction ([9T|) . 

1 1 z „ / rxs \ rxs 

6 9z = ^hr U(X5) + e r | Xs -v s -2 / dx* / d X 2 * . (100) 
X s H s V Jo J Jo 

The observed reduced shear at constant redshift becomes I be 

2 V 2 £> 5V 

* = -2jJ + (W +% ' <101) 
where the first, the second and the third terms on the right-hand side are respectively given by 
eqs. ([72]) . ([95]) and ([lOOp . This is the main result of this article. 



8 Note that the only effect of the transverse velocity of the source is to modify the direction under which 
the galaxy is observed. Thus, it has no effect on the shear. 
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5 Conclusion 



In this article we have derived the expression of the reduced cosmic shear up to second order in the 
perturbations with full-sky validity. Our main result is summarized in eq. (jlOip . As it is expressed 
in terms of spin-2 operators on the sphere it can be decomposed as sum of spin-weighted spherical 
harmonics on the sky. Indeed, this description ensures that our observable has a genuine spin-2 
behavior on the celestial sphere. 

Our result is written in terms of the metric perturbations in the generalized Poisson gauge. These 
are the scalar potentials <j> and ifi and the vector and tensor components of the metric generated at 
second order, respectively ui{ and hij. Let us first comment on the first two terms on the right-hand 
side of eq. ()101|) . Remarkably, the contribution from scalar perturbations from the sum of these two 
terms can be expressed in terms of the Weyl potential \& = (0 + ^)/2 only. As explained, this is due 
to the fact that null geodesies are conformally invariant. These two terms contain the well-known 
second-order corrections due to lens-lens coupling and departure from the Born approximation, 
which dominate in the small-angle approximation. On larger angular scales new couplings become 
important. These are an intrinsic contribution which is a purely general relativistic effect at second 
order, a coupling between the gravitational potential at the source with the lens and corrections due 
to couplings between the lens and the photon time-delay. We have checked that these contributions 
can be independently reconstructed from the calculation of the shear at first order in a universe with 
a radially dependent spatial curvature. Other checks, such as the invariance under a homogeneous 
time shift and a conformal transformation can be used to verify the validity of these new corrections. 
Another scalar correction appears in the form of products of two spin-1 fields and comes from the 
couplings between two photon deflections. Finally, besides the scalar contributions, the shear gets 
a contribution from spin-2 quantities defined from the vector and tensor components of the metric 
generated at second-order. Note that the separation between all these contributions is not gauge 
invariant. 

In Poisson gauge, the correction due to the coupling between the photon redshift perturbation 
and the lens cannot be written in terms of Vl/ only. Indeed, the integrated contribution to the photon 
redshift - the integrated Sachs- Wolfe effect - is a time integral over ^ but the intrinsic contributions 
- Sachs- Wolfe and Doppler effects - are expressed in terms of the Newtonian gravitational potential 
and the velocity along the line of sight and do not depend on ^/ only. 

We are now in the position to explore the phenomenological consequences of these results in view 
of the future (partially) full-sky lensing surveys. In particular, the new corrections that we have 
computed should become relevant in deriving the lensing bispectrum on large angular scales. For 
instance, to compute the bispectrum in the squeezed limit one needs to take one of the three modes 
to be much smaller than the other two, corresponding to angular scales comparable to the depth 
of the survey. As the lensing is a cumulative effect integrated along the line of sight, it is difficult, 
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at this stage, to precisely guess the relative importance of the various contributions. In particular, 
although the lens-lens coupling terms are a priori larger by a factor ~ I 2 compared to the others, 
they may be damped by geometrical factors. We leave these investigations for the future. 
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A Conformal relations 

We consider two conformally related metrics g a p = ^l 2 (x^)g a p. We want to show that the cor- 
responding Jacobi mappings computed using eq. (fTUj) with the respective metrics are related by 
D ab = T> ab . We can decompose the Riemann tensor in terms of the Ricci tensor R a p and the Weyl 
tensor C a p pa as [34] 

where the brackets in the indices denote antisymmetrization. With this decomposition, the definition 
of TZ a b, eq. Q, simply yields 

Tl ab = C a ^ 5 n a a kPtfn b s - ~5 ab R^kPk"< , (103) 

where we have used the normalization properties of nj 1 and k^. 

Using this expression, let us study how lZ a b transforms under a conformal transformation. As 
ria is normalized to unity it transforms as ha = Q~ 1 nJ 1 . If A is an affine parameter of the null 
geodesic in the metric g a p, then the affine parameter computed using the metric g a p is related to 
A by dA = Jl 2 dA [34 . Thus, = S1 -2 A^. Furthermore, as the Weyl tensor with one upper index 
Ca/35 7 is invariant under conformal transformation, the first term of (|103p transforms as 

C a ^ s h a a ¥frn b s = n- A C aMS n a a k^k^n h s . (104) 

The Ricci tensor transforms as [33] 

R ai = R ai -2V a V 1 lnn-g ai g Sa V s V a lnn + 2(V a ln{l)V J lnn-2g aj g 5a (Vs^n)V a lnn . (105) 
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Projecting this expression by k a k' ) , the two terms proportional to g ai vanish because of the null 
condition of the photon wave vector, while the covariant derivatives can be written as derivatives 
along the null geodesic. Using the geodesic equation one obtains 

2" 



- 2 



d 2 lnfi 



+ 2 



din SI 
dA 



Thus, from eq. ()103j) lZ a b transforms as 



d 2 lnil 
^ — dA 2- 



dlnft 
dA 



Finally, using this transformation it is easy to show that if V a b is the solution of 

d 2 



dA 2 



V. 



ab 



KarV, 



ac^cb i 



then V a b = Q, V a f, is the solution of the corresponding equation for the metric g^, 

dA 2 



(106) 



(107) 



(108) 



(109) 



The relation between and T> a \, can be easily understood by noting that V a i, relates distances 
at the source, that scale like 0, to angles at the observer that are invariant under a conformal 
transformation. 



B Spin operators on the sphere 

The construction of spin fields can be easily done on a plane identified with the complex plane of 
coordinates z = x + iy. Let us consider a complex field s X(z) whose value depends on z. This 
field will be said to have spin-s if its value is changed in e lsa s X(z) after a rotation of angle a. For 
instance, if qX{z) is a scalar (i.e. a spin-0) field, then 

x x(z) = d x x(z) + id y x(z) , (no) 

is a spin-1 field. This relation can alternatively be written as \X(z) = 2d z qX(z), where the partial 
derivative is to be taken for a fixed value of z. In general, the successive application of the operator 
= 2d z leads to the construction of spin-s fields. Equivalently, the operator = 2dz lowers the spin 
by one. In the context of standard first-order lensing theory, the complex shear field 7 = 71 + i72 is 
a spin-2 field that derives from the projected potential ip, i.e. 7 = 2 ip. 

This construction can be extended to the sphere when one does not want to rely on the small- 
angle approximation. The early elements of such a construction date back to [42|, 140] . In general, one 
is naturally led to introduce the Euler angles (9, ip, a) so that the coordinates of a point on the unit 
sphere are (sin# cos<^, sin 9 sin ip, cos 9). To each of these points one can associate a radial vector e r 
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and two tangential vectors ei and &2 that can be conveniently chosen along the lines p = const and 
9 = const respectively if a = 0. The angle a then corresponds to a rotation around the axis e r that 
rotates eg and e v with an angle a. 

As for the plane, a spin-s field is such that its phase varies as sa, by rotation of an angle a. The 
construction of the operators and relies on the use of the complex stereographic coordinates, 



which map the sphere onto the complex plane. More specifically it can be shown that the expressions 
of the operators and depend explicitly on the spin of the field S X to which they are applied, 



S X = -(sinfl) 5 {d e + icsc0d v ) (sm9)-%X and S X = -(shitf)" 5 (d e - icscfld^) (sm9) s s X . 

(113) 



In analogy to the case of the plane, the Laplacian operator formally reads A = 0. Note however 
that the relation 00= 00 holds only when the operators act on a scalar field. In general we have, 
(00- 00) s X = 2s s X. 

Spherical harmonics, Y" z m (#,</3), are spin-0 functions that are the eigenfunctions of the Laplacian 
with eigenvalue -1(1 + 1), i.e. AY l m (9,p) = -/(/ + l)Y l m (9,p) with 8^(9, (p) = imY, m (0,<^) and 
with a specific normalization. The orthogonality relation, 



is the key property that makes it possible to decompose any function into spherical harmonics. 

In general spin-s fields are decomposed on the basis of the spin- weighted spherical harmonics, 
that can be obtained through the application of the operator and on the spherical harmonics. 
More specifically, we define s Y l m (9, ip) with 




(111) 



S X = 2P 1 - 3 d c P s s X and S X = 2P l+s d^P~ s s X , 
where P = ^(1 + CO- Expressing these operators in terms of 9 and p, one finds 



(112) 




(114) 




(115) 



(116) 



The spin- weighted spherical harmonics obey the following relations, 



s Yr(o,p) 

00 s Yr(e,p) 



-(i-s)(i + s + i) s Yr(e,p) 



(117) 
(118) 
(119) 




(120) 
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with J d 2 Q = Jq W dip J_ 1 dcos 6. 

One usually defines two scalars, E and B, associated to the spin-2 shear 7. To do that, one 
decomposes the shear 7 and its complex conjugate 7* as, 

1 (6,<p) = ^2 2 a lm2 Y l m (6,<p) and j*(6, <p) = £ _ 2 a im _ 2 ^ m (0, <p) . (121) 

im im 

As in the context of CMB polarization (see [SI]), E and B can be defined through their harmonic 
decomposition, 

E(0, <p) = ~^2( 2«im + - 2 ai m ) Yr(9, ip) and B(6, tp) = i £ ( 2 a im - _ 2 a im ) ^ m (0, ^). 

im im 

(122) 

I? is invariant under parity change, whereas B changes signs. 
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